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O ■ Large- A'^ expansions and computer simulations indicate that the universality class of 

Q\ • the finite temperature chiral symmetry restoration transition in the 3D Gross- Neveu model 



j^ I is mean field theory. This is a counterexample to the standard 'sigma model' scenario which 

predicts the 2D Ising model universality class. We trace the breakdown of the standard 
scenario (dimensional reduction and universality) to the absence of canonical scalar fields 
in the model. We point out that our results could be generic for theories with dynamical 
^ ■ symmetry breaking, such as Quantum Chromodynamics. 
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When studying the finite temperature chiral restoration transition in QCD one is 
usuaUy guided by the concepts of dimensional reduction and universahty. A compeUing 
idea, first put forward in [1] and later elaborated in [2], is that in four-dimensional QCD 
with Nf light quarks the physics near the chiral transition can be described by the three- 
dimensional cr-model with the same global symmetry. The reasoning behind this proposal 
is based on counting the light degrees of freedom and can be phrased as follows. The 
transition region is dominated by the longitudinal and transverse fiuctuations of the order 
parameter, a and tt, which go soft at the transition temperature. Being bosonic, a and tt 
have zero modes, uj^ = 0, in their finite-temperature Matsubara decomposition. These zero 
modes are the only relevant degrees of freedom in the scaling region and at low energies the 
n 7^ modes decouple. Therefore, in the context of a (i-dimensional theory, one concludes 
that the phase transition is described by an effective scalar theory in d—1 dimensions. As 
a consequence, the chiral transition of four-dimensional QCD, with Nf = 2 fiavors, should 
lie in the same universality class as a three-dimensional 0(4) magnet [1,2]. Similarly, 
other models e.g. four-fermi theories in d-dimensions like Gross-Neveu [3] with discrete 
or Nambu-Jona-Lasinio [4] with continuous chiral symmetries, are expected to be in the 
universality class of a d — 1-dimensional Ising or Heisenberg magnet, respectively. 

It is the purpose of this paper to discuss the assumptions underlying this analysis. As 
an illustration we will study two examples: a purely bosonic theory, 0{N) a-model where 
the ideas of dimensional reduction apply, and a Gross-Neveu model with composite scalars 
where they fail. We discuss the generic features of the models that might apply to other 
field theories at finite temperature. At the end we comment on the implications these two 
examples have on QCD. 

To illustrate how the idea of dimensional reduction is realized in scalar theories, we 
start with the A^-component scalar theory and consider the large- A^ limit [5] for simplicity. 
To avoid complications due to Goldstone bosons, we work in the symmetric phase. At zero 
temperature, the susceptibility is given by the single tadpole contribution. Defining the 
critical curvature fx'^ as the point where the susceptibility diverges (/U^ + -^/ 1/?^ = 0)? 
the expression for the inverse susceptibility can be recast into 
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where / = J d'^q/ {2tvY, and we absorb the combinatorial factor in A. The extraction of 
the critical index 7 reduces to counting powers of the infrared (IR) singularities on the 




left hand side (LHS) of eq.(l). Above four dimensions, both terms are IR finite and the 
scahng is mean field (7 = 1). Below four dimensions the second term in eq.(l) dominates 
the scaling region - the integral diverges as ■)(^^~'^>''^ . This gives the zero-temperature 
susceptibility exponent 7 = 2/((i — 2) [5]. 

At finite temperature, apart from the replacement of the frequency integral with the 
Matsubara sum, modifications are minimal [6] . For a given value of /U^ we define the critical 
temperature, Tc, by /U^ + ATc Xln /- l/l^^nc + ^) = 0, where Unc = 2TvnTc. The momentum 
integrals are now performed over d — 1 dimensional space. Separating the n = mode 
{ujq = 0) from the rest of the sum, we get the leading singular behavior 

The n = piece dominates the scaling region. It resembles the zero-temperature expres- 
sion, eq.(l), except that now, the integrals are performed in (i — 1 dimensions, instead of d. 
The power counting is the same as before and it yields the thermal exponent 7^ = 2/((i— 3) 
which is the same as the zero-temperature 7 in (i — 1 dimensions [6] . It is easy to obtain 
the other critical exponents; they show the same type of behavior as 7. 

To illustrate how compositeness affects the physics near the phase transition, we an- 
alyze the problem of chiral symmetry restoration in a Gross-Neveu model given by the 
lagrangian L = ijj(id + m + gcr)tjj — |(t^, where notation is standard [3]. Besides being an 
interesting theoretical model, it is also believed that, when properly extended to incorpo- 
rate continuous chiral symmetry, four-fermi models are more realistic as effective theories 
of QCD than the linear sigma model, especially at scales where quark substructure is im- 
portant. When fermions are integrated out of the Gross-Neveu model, the Ising symmetry, 
a -^ — (J, of the effective action becomes manifest. If the dimensional reduction -t- uni- 
versality arguments hold [1] , the finite temperature transition of the (/-dimensional model 
would lie in the universality class of the d — 1 dimensional Ising model. In the remainder 
of the paper we explain how and why this argument fails. 

First, we start with the zero-temperature gap equation and corresponding critical 
exponents. The model can be treated in the large- A^ limit. To leading order, the fermion 
self-energy, E, comes from the cr-tadpole: E = m — (7^ < t/jtl; >. To obtain the scaling 
properties of the theory, we define the critical coupling as 1 = 4g'^ J l/q^. Combining this 
definition with the gap equation leads to 
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Like the scalar example, this form is especially well suited for extracting critical indices 
since the problem reduces again to the counting of the infra-red divergences on the right 
hand side [7,8]. The critical indices are defined by < ipijj > \m=o ~ ^^7 < V'V' > |t=o ~ 
m^/^,E|^=o ~ t^i etc.. Here, t = g'^/gc — 1 is the deviation from the critical coupling. 
Since S ~< t(jt(j >, /3 = z/ to leading order. Above four dimensions the integral in eq.(3) is 
finite in the limit of vanishing E and the scaling is mean-field. 

Below four dimensions, the E — > limit is singular - the integral scales as E*^"^. Thus, 
in the chiral limit t ~ E'^"^, and at the critical point, t = 0, away from the chiral limit, 
m ~ E'^"^. The resulting exponents are non-gaussian: (3 = l/{d — 2) and d = d — 1. The 
remaining exponents are obtained easily: ?7 = 4 — (i, 7= 1 [7,8] and one can check that 
they obey hyperscaling. 

We now consider the Gross-Neveu model at finite-temperature. We choose to stay 
between two and four dimensions to emphasize how zero-temperature power-law scaling 
changes at finite temperature. The gap equation is now modified to 

where Un — {2n + l)nT. For g > gc the critical temperature is determined by: 1 = 
4Tc(7^ Xln /^-'-/('-^nc + 9^)5 where iVnc = (2n + 1)ttTc. This expression defines a critical 
line in the {g, T) plane. For every coupling there exists a critical temperature beyond 
which the symmetry is restored. Conversely, for a fixed temperature there is a critical 
coupling, defined by the above expression, corresponding to symmetry restoration. At 
zero temperature, the symmetry is restored aX g = gc- Thus, {g = g^T = 0) is the 
ultra-violet (UV) fixed point. As the coupling moves away from g^ a higher restoration 
temperature results. At infinite coupling the end-point, ((7 = oo,T = Tc), is the IR fixed 
point. The critical line connects the UV and IR fixed points dividing the {g,T) plane into 
two parts. ^ 



^ The equation for the critical line can be brought into a compact form by combining 
the expression for T^ with the definition of the zero-temperature critical coupling. This 
results in: [g"^ / g1 — 1) ~ T^~'^{g), i.e. Tc{g) ~ E(T = 0). In this way, for any value of the 
coupling, the critical temperature remains the same in physical units. 



Combining the definition of Tc witfi tlie finite-temperature gap equation, we can bring 
it to a form similar to eq.(3) 
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Tlie extraction of tfie critical exponents proceeds along the same lines as in the zero- 
temperature case. One difference relative to eq.(3) becomes apparent immediately: the 
zero modes are absent here and the integrand in eq.(5) is regular in the E ^ limit even 
below four dimensions. Consequently, the IR divergences are absent from all the integrals 
and the scaling properties are those of mean-field theory: P = u = 1/2, 6 = 3, etc. This 
is true for any d, below or above four. It appears that in this case, contrary to the scalar 
example, the effect of making the temporal direction finite (1/T) is to regulate the IR 
behavior and suppress fiuctuations. This is manifest in other thermodynamic quantities 
as well. For example, to leading order, the scalar susceptibility, x = 9 < ?/'?/; > /dm, is 
given by 

Once again, because of the absence of the zero mode (cuq = t^T), the integral in eq.(6) 
is analytic in E, and the mean field relation x~^ ~ ^^ follows. This is equaivalent to 
7 = 2z/ = 1. The explicit calculation of the momentum dependence of the a propagator 
[9] yields ?] = 0. 

The scaling laws of the finite temperature transition obtained above are completely 
different from the predictions of ref.[l]. In fact, even the systematics are opposite. The 
regulating character of the temperature drives the lower dimensional theory towards an 
effective theory that has gaussian critical exponents. 

The fermionic model discussed above was first analyzed in ref. [9] . Higher order calcu- 
lations have shown that the results are not artifacts of the large- A^ limit [10]. In addition, 
it was explained in [10] how the Ising point is recovered in four dimensional Yukawa models 
beyond the leading order in 1/A^ and why this does not happen in Gross Neveu models. 
Lattice simulations of the three dimensional model have verified the predictions of the 
large- iV expansion at zero temperature, at nonzero temperature and at nonzero chemical 
potential [11]. The results for critical indices have been verified and improved by larger 
scale simulations enhanced by histogram methods [12]. We have done additional simu- 



lations to check the finite temperature results [11] in detail. Lattices of sizes 6 x 30^, 
12 X 36^ and 12 x 72^ were simulated at A^ = 12 using the Hybrid Monte Carlo algorithm 
described in [11]. High statistics runs (several tens of thousands of trajectories for each 
coupling) were made on a variety of lattices to guarantee that the simulations were probing 
the physical IR modes at finite temperature. Luckily, our task is to distinguish mean field 
exponents from those of the two dimensional Ising model, and, as reviewed in Table I, 
they are dramatically different. We will discuss the exponents d and (3 (defined above) 
and leave other calculations to a lengthier presentation. In Fig.l we show the square of 
the order parameter a plotted against 1/g'^. The data is in excellent agreement with mean 
field theory where P = 1/2, and rules out the Ising model value of 1/8. Note that the 
statistical error bars in the figure are smaller than the plotting symbols themselves. Since 
both lattice sizes give the same estimates of /3 while their critical temperatures are quite 
different, we are confident that the simulation is probing the true continuum behavior 
of the finite temperature transition and is not corrupted by a sluggish crossover between 
symmetric and asymmetric lattices. Several runs on the huge 12 x 72^ lattice were made 
and the 12 x 36^ results were confirmed to a fraction of a percent. We also calculated the 
susceptibility index 7 and found 7 = 1.0(1) in the same runs. The Ising result 7 = 7/4 
is decisively ruled out. Next, we read off the critical temperatures for both lattice sizes, 
measure the response of the order parameter at criticality to an external symmetry break- 
ing field ( bare fermion mass ) and obtain 6. The data is shown in Fig. 2, and for both 
lattice sizes we find 6 = 3.1(1). The Ising model value of 5 = 15 is ruled out. In all of these 
calculations we carefully visualized the a field to check for nonuniform configurations that 
would violate the mean field hypothesis [13]. None were found and all the past simulations 
[11] and the new ones reported here support the contention that the large- A^ results are 
reliable for this problem. 

An important feature of the exponents corresponding to the finite-temperature tran- 
sition, is that they violate hyperscaling [14]. Usually, hyperscaling violation occurs above 
four dimension and is expressed in terms of exponent inequalities [14] e.g. 2/35 — 7 < du. 
Strict inequality is applicable only for d > 4 and implies factorization of the correlation 
functions. In our example the above inequality goes in the opposite direction and the 
breakdown of hyperscaling is not accompanied by the factorization of Green's functions. 

In conclusion, the study of the Gross-Neveu model suggests that arguments invoking 
dimensional reduction + universality must be used with care. Our results indicate that 
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an effective scalar model fails to describe the Gross-Neveu model at finite temperature. 
We believe that the reason for this failure in our example is related to the composite 
nature of the mesons. Pointlike scalars can not adequately describe the physics in the 
vicinity of the second order chiral transition. The physical picture behind this failure 
observes that both the density and the size of the loosely bound sigma meson increase with 
temperature. Close to the restoration temperature the system is densely populated with 
overlapping composites. In other words the fluffiness of the mesons can not be ignored 
- the constituent fermions are essential degrees of freedom even in the scaling region, 
right before the composites dissociate. Similar discussions of the failure of effective meson 
theories in a slightly different context have been given in ref.[15]. 

It is well known that four-fermi models can be used as effective theories of QCD 
[16]. In addition to having the same global symmetries, the mesons in both theories are 
composite. Therefore, these models are believed to have common properties over a wide 
range of scales where the quark substructure of the mesons is relevant. Knowing this, 
it would be interesting to see what happens in two flavor QCD [17]: does it follow the 
dimensional reduction scenario, or the Gross-Neveu behavior? Of course, these two alter- 
natives do not exhaust all the possibilities [18], but we believe that the scenario suggested 
by the Gross-Neveu model is sufficiently compelling to warrant further analyses of QCD 
simulation data. 

We wish to acknowledge the discussions with Eduardo Fradkin and Maria-Paola Lom- 
bardo. Special thanks go to Rob Pisarski for his comments on the manuscript. This work 
is supported by NSF-PHY 92-00148 and used the computing facilities of PSC and NERSC. 



References 

[I] R. Pisarski and F. Wilczek, Phys. Rev. D29, 338 (1984). 

[2] F. Wilczek, Int. J. Mod. Phys. A7, 3911 (1992); K. Rajagopal and F. Wilczek, Nucl. 

Phys. B404, 577 (1993) 

[3] D. Gross and A. Neveu, Phys. Rev. D20, 3235 (1974). 

[4] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 

[5] See, for example, C. Itzykson and J.-M. Drouffe, Statistical Field Theory (Cambridge 

University Press, 1989). 

[6] L. Dolan and R. Jackiw, Phys. Rev. D9, 3320 (1974). 

[7] J. Zinn- Justin, Nucl. Phys. B367, 105 (1991). 

[8] S. Hands, A. Kocic and J. B. Kogut, Phys. Lett. B273 (1991) 111. 

[9] B. Rosenstein, B. Warr and S. Park, Phys. Rev. D39, 3088 (1989). 

[10] B. Rosenstein, A. Speliotopoulos and H. Yu, Phys. Rev. D49, 6822 (1994). 

[II] S. Hands, A. Kocic and J. B. Kogut, Ann. Phys. 224, 29 (1993); Nucl. Phys. B390 
355 (1993). 

[12] L. Karkkainen, R. Lacaze, P. Lacock, and B. Petersson, Nucl. Phys. B415[FS], 781 

(1994). 

[13] R. Dashen, S. K. Ma and R. Rajaraman, Phys. Rev. Dll, 1499 (1975); F. Karsch, J. 

Kogut and H. W. Wyld, Nucl. Phys. B280[FS18], 289 (1987). 

[14] B. Freedman and G. A. Baker Jr, J. Phys. A15 (1982) L715. 

[15] R. Jaffe and P. Mende, Nucl. Phys. B369, 189 (1992). 

[16] See for example a recent review by T. Hatsuda nad T. Kunihiro, Tsukuba preprint, 

UTHEP-270, (to appear in Phys. Rep.). 

[17] F. Karsch, Nucl. Phys. (Proc. Suppl.) B34, 63 (1994). 

[18] E. Shuryak, Comments Nucl. Part. Phys. 21, 235 (1994). 

Figure Captions 

1. Order parameter squared plotted against temperature on 6x30^ (left) and 12x36^ (right) 
lattices. 

2. Order parameter response at criticality plotted against bare fermion mass on 6 x 
30^ (bottom) and 12 x 36^ (top) lattices. 
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Critical exponents of the 3D Gross-Neveu and 2D Ising model 
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Order Parameter Response at Criticality, 6x30^2 and 12x36^2 
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